Whence the odd-even staggering in nuclear binding? 
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We explore the systematics of odd-even mass staggering with a view to identifying the physical 
mechanisms responsible. The BCS pairing and mean field contributions have A- and number parity 
dependencies which can help disentangle the different contributions. This motivates the two-term 
parametrization ci -|-C2/^ as a theoretically based alternative to the inverse power form traditionally 
used to fit odd-even mass differences. Assuming that the A-dependence of the BCS pairing is weak, 
we find that mean-field contributions are dominant below mass number A ~ 40 while BCS pairing 
dominates in heavier nuclei. 



I. INTRODUCTION 

The observed odd-even mass staggering has been ex- 
tensively explored in the literature The eflFect 
has generally been associated with the pairing gap A, as 
suggested by BCS theory. When examined over the full 
table of measured nuclides , the values for the stagger- 
ing display a great deal of scatter. However, a slow trend 
with nuclear mass number A is also observed. This trend 
has been extensively analyzed in the literature by fitting 
data to an inverse power of A. Traditionally this has 
been 1/A^^^, but other powers are used as wellfl, 0, [1]. 
For example, ref. [3| advocates an expansion in powers of 
A"^/^, based on the finite-size analysis in ref. [i7|- There 
is no theoretical justification for a single-term fractional- 
power dependence, however. 

It has been suggested that, in addition to the antici- 
pated BCS type of pairing, 

other effects contribute significantly to the staggering 
P) H, IE]- One would anticipate that each contribution 
would have its own characteristic dependence on A. Con- 
sequently the staggering should follow some combination 
of these. The tradition of representing the net effect by a 
form with a single fractional exponent obscures the com- 
plexity of the effect. We propose here a different form 
that may offer a better opportunity to identify the vari- 
ous contributing effects. 

There are a number of difference formulas that may 
be used as measures for the staggering 0, H, In this 
work we use the measure conventionally called A^'^^ [H, 0, 
01 • It involves the difference between the two separation 
energies of adjacent nuclei. Here the removal of one type 
of nucleon, i.e., proton or neutron, is considered, while 
the number of the opposite isospin is kept the same. We 
believe it is the most direct measure of staggering and is 
the most convenient quantity for theoretical study. One 
defines the neutron OES as 

A(3)(iv, z)) = -i-r- i {Sn{N + 1,Z)- Sn{N, Z)) 

(1) 

Here Sn is the separation energy, Sn{N) — B{N)—B{N— 
1) and B the binding energy. The factor depending on 
the number parity tt^t is chosen so that the OES cen- 



tered on even and odd neutron numbers N will both be 
positive for normal attractive pairing. The proton OES 
are defined similarly and we designate them by Aop and 

A(3) 

The 4-point measure used by Bohr and Mottelson 
Eqs. 2-92, 2-93] is the average of adjacent even and odd 
values of A^^^^ . We will also find it convenient to examine 
the average and differences of our parametrized OES, but 
that does not require contiguous data. It should also 
be noted that, because the A^'^) values represent second 
differences of the binding energies, they may be non-zero 
even without any odd-even staggering. 

The various OES offer interesting possibilities for ex- 
amining the various contributions. First, the even and 
odd OES have different contributions that can possibly 
be distinguished. Second, the two isospins (neutron and 
proton) can be compared to give some indication of the 
role of the Coulomb and other isospin dependent effects. 



II. MECHANISMS 

There are a several mechanisms that contribute to OES 
in framework of mean field theory. We consider three of 
them: single-particle energies, BCS pairing, and orbital 
interaction energies. We will assume that orbitals are 
two-fold degenerate between time-reversed partner states 
to estimate single-particle and diagonal interaction ener- 
gies. A summary of the following discussion is provided 
in Table I below, which lists separately the contributions 

(3) (3) 

pertinent to Ao and Ae . There are several kinds of or- 
bital interaction that affect the systematics of the OES. 
Most important are the spin- and symmetry-dependent 
interactions that are classified under "time-odd fields" 
in the density functional formulation of mean-field the- 
ory. But also the Coulomb interaction and the ordinary 
isospin-spin dependent interactions play a role when ex- 
amining the differences of A^^") measures. Another mech- 
anism that has been suggested [U, 0] involves the static 
polarization of the core by the valence orbitals. However, 
it has no obvious parametric dependence on A, and we 
must assume that its contribution is small to carry out 
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the present analysis. 



TABLE I: Contributions to odd-even staggering of nuclear 
binding energies, and their assumed yl-dependence. 



Mechanism 




Af' 


A-dependence 


Single-particle 





(ei - ei-i)/2 


l/A 


BCS correlation 




^BCS 


const. 


time-odd 






l/A 



A. Single-particle energies 

This mechanism is discussed in many places, e.g., ref. 
p'l, and we repeat the argument for completeness. The 
binding energy associated with the single-particle Hamil- 
tonian is the negative of the eigenvalue e^. The neutron 
(or proton) orbitals are doubly degenerate implying that 
the contribution from the two separation energies odd 
OES. On the other hand, the two separation energies in 
the even OES come from different orbitals. The resulting 
OES are 

A(3)(s.p.)=0; A(3)(s.p.) = (e.-e._i)/2. (2) 

For brevity we will refer to this as the "single-particle" 
contribution. The average difference in single-particle 
energies can be estimated in the Fermi gas model as [9^ 



^ep/iN « 93/A, MeV. 



(3) 



In the last approximate equality we have estimated the 
level spacing using Ci? « 35 MeV and A « 2N . 



B. BCS correlation energy 

To avoid confusing the BCS correlation energy with 
the mean-field energies, let us define it as the additional 
energy associated with an interaction of the form 



(4) 



The orbital i is the time-reversed partner of orbital i. The 
prime on the summation indicates that a sum over i takes 
only one member of the partners («,«). The restriction to 
i ^ j ensures that there is no mean-field contribution. In 
the large number limit, the BCS energies give the familiar 

result, Ao,e = /S.Bcs^ where Abcs is the pairing conden- 
sate energy for an orbital at the Fermi surface. We will 
refer to this contribution as the "pairing" contribution. 

The A-dependence of the BCS pairing is an open ques- 
tion at present. An analytic treatment can be carried 
out[l3| following the derivation of the BCS formula, 



Here S is the half-width of the truncation zone, G is an 
average pairing interaction matrix element and is 
the density of levels at the Fermi-energy, pi. For large 
A, we have G ^ l/A and p ^ A giving a constant A. 
At the next order in a finite-size expansion both G and 
p add terms varying as A~^/^. This give s an overall A- 
dependence A{BCS) ~ exp{c/A^/^) However, the 

finite range of the nuclear interaction as well as the in- 
duced pairing introduces additional A dependencies that 
are hard to estimate. In a recent global study fl§\ it 
was found that the A dependence of the average OES is 
quite weak, taking pairing interactions with density de- 
pendence as commonly used. We shall therefore make 
the simplest assumption here, that the BCS correlation 
energy contribution is independent of A. Finally, we note 
that although Eq. ^ was derived only for the Hartree- 
Fock limit, it remains approximately valid in the presence 
of the pairing Hamiltonian Eq. 



C. Time-odd fields 

A semi-analytic treatment of this mechanism is given 
ll| . There it is called a polarization effect, but in 



ref. 



fact it is present in the diagonal orbital interaction. The 
expressions for the effect are quite different depending on 
whether one is using a density functional formalism or a 
formalism based on interaction matrix elements between 
orbitals. The first formalism has unphysical self-energies 
that are canceled by the exchange term in the interaction. 
In the second formalism self-interactions do not appear 
and one does not make use of single-particle fields and 
their transformation properties under time reversal. For 
simplicity in the following we follow the approach based 
on interaction matrix elements. In this formulation, we 
count the interactions as particles are added to the sys- 
tem, always filling both orbitals of a Kramers pair if pos- 
sible. Suppose we start with an even-even system with 
A nucleons and add another neutron in orbital i. The 
additional interaction energy is 



IE' 



Av. 



Here the j summation runs over all occupied orbitals ex- 
cept i and we have defined an average interaction energy 
V in the last equality. We next add another particle to 
the system in orbital i, giving an additional interaction 
energy, Vi + Vfi n- Then Ao'^-' is given by half the differ- 
ence, 



A(3) = 



^/2 



(6) 



A{BCS) = 2S'e-i/(<=^''(^». 



(5) 



Adding a third particle gives an additional interaction 
energy of {A + 2)v. This can be used to obtain the fol- 
lowing estimate for Ae 

Af^=v- v^^ii/2. (7) 
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For the particle-particle interaction matrix elements as- 
sociated with the time-odd contribution we use the same 
argument as with the pairing matrix elements G. This 
leads to the expectation that the average matrix elements 
embodied in the energies v and ^ also would have an 
dependence. The formulas Eq. ([6]) and (I7| for the 
time-odd contribution were derived assuming that the 
orbital occupancies are zero or one, but the formulas re- 
main valid in the presence of pairing in the nonblocked 
orbitals, as was the case for the single-particle contribu- 
tion. 



TABLE IL Fits to the OES. Column 2 reports the size of each 
data set. Columns 3-5 show the rms residuals for the three 
fit functions discussed in the text. Columns 6-8 gives the fit 
parameters for models i) and ii). Units are MeV. 



Data set 


rms residual 


fit parameters 


OES 


Size 


i ii Hi 


Co 


Cl C2 


^(3) 


443 


0.313 0.254 0.270 


1.04 


0.82 24. 


^(3) 


442 


0.420 0.278 0.308 


1.32 


0.94 41. 


^(3) 


418 


0.275 0.220 0.231 


0.96 


0.75 25. 


^(3) 


407 


0.455 0.245 0.270 


1.64 


1.11 59. 



III. EXAMINING THE EXPERIMENTAL DATA 

We now examine the experimental data on the OES, 
fitting it to various parametrizations. The data sets are 
defined with the same criteria used in ref. p^ . using 
the Audi 2003 mass table 0. We only include nuclei for 
which the experimental error on the binding energy is 
less than 200 keV, and make a further selection on the 
experimental data to avoid special binding effects not 
related to the mechanisms we have discussed. We drop 
nuclei with N = Z because of their Wigner energy jlj, 
[T^ contribution. We also drop odd-odd nuclei which have 
an additional neutron-proton pairing p^ . [l^ . Finally, to 
isolate isospin dependencies, we require that the nuclei 
have positive isospin {N > Z) . The sizes of the resulting 
data sets are given in Table II. 

For each of the data sets we fit the OES with the fol- 
lowing fitting functions: 

i) A constant cq. This is the naive BCS form, and sets 
the scale for any improvements of the phenomenol- 
ogy. 

ii) The one-term, two-parameter function c^A^"' with 
prefactor Cq and exponent a which provides the 
minimum rms deviation. This shows the best one 
can do with two parameters, irrespective of justifi- 
cation. 

iii) The two-term expression ci + cij A. The constant 
Cl represents the BCS correlation energy, assumed 
independent of A, and the term CilA~^ represents 
interaction and single-particle contributions. We 
use additional subscripts o, e and n,p to distinguish 
the different OES measures. 

The rms residuals of each of the four OES's with each 
of the three fit functions are shown in Table II. We also 
show the fit parameters for fits i) and iii). 

A. Quality of fits 

Let us first examine the quality of the fits. The sim- 
plest, just a constant, has an rms residual of about 0.3 
MeV for the odd OES and 0.45 for the even OES mea- 
sure. The residuals are decreased substantially, by up 



to 46%, going to the power law fit, model ii). Model 
iii) with two terms decreases the rms residual up to 41% 
and are almost competitive with model ii) . In view of its 
better theoretical justification, we would advocate using 
it instead of the power-law fit. We note that the rms 
residual is larger for the even OES than for the odd one. 
This already suggests that shell fluctuations associated 
with the single-particle mechanism may be stronger than 
those affecting the time-odd contribution. 



B. Neutron fit parameters 

We now examine what information can be extracted 
from the fit parameters that we have obtained in model 
iii). Examining the neutron OES, we first note that 
the constant terms in the even and odd OES are close 
to each other; the ratio of the difference to the sum 
is 0.06/0.88 « 7%. Since the BCS is the only A- 
independent contribution, and it is common to both the 
even and the odd OES, the near-equality is just what we 
expect. Comparing the values to the one-term constant 
Co ~ 1.0 — 1.3, we see that the assigned BCS pairing 
represents only ~ 75% of the total on average. 

The coefficients ci , which we associate with the mean- 
field mechanisms, are quite different in the even and odd 
OES. That is also anticipated, since the two mechanisms 
we discussed affect them differently. From the coefficients 
C2o,C2e and an assumed single-particle contribution ac- 
cording to Eq. ([3]) , we can extract average values for ma- 
trix elements in the time-odd contribution. These values 
are: 

WiJ.ij ~ -49 MeV; u«-30MeV (8) 

The signs and relative magnitudes are as expected for a 
basic interaction is attractive and short-ranged. For 
such interactions, the exchange contribution is large and 
attractive for V{{ {i but not for v. However, we cannot ex- 
pect to extract reliable values from this simple approach 
for reasons already mentioned. 

It is useful at this point to average the even and odd 
OES to emphasize the BCS correlation contribution and 
to look at the differences as well which will empha- 
size the mean field contribution. We define parameters 
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Ck± = {ckoicke)/'2 and show the fits in this form in Table 
III. The entries for the averages show again the relative 
importance of the mechanisms. The two terms give equal 
contribution at mass A satisfying Ci+ = C2+/A, giving 
A ~ 40 as the dividing point between light nuclei domi- 
nated by non-BCS pairing and heavier nuclei dominated 
by BCS pairing. We have already remarked that Ci_ is 
small, as expected from a BCS pairing mechanism. The 
Cfe_ arc essentially determined by the curvature of the 
binding energy surface in the two directions, along the 
Z axis for c^-p and the N axis for Cfe_„. In the liquid 
drop model, this curvature is mostly due to the Coulomb 
and symmetry terms in the formula. The curvature can 
be largely fit with the C2-, as expected from the func- 
tional form of the symmetry term, ~ (A^ — Z)'^/A. The 
Coulomb is different, however, as we discuss below. 

TABLE III: Averages and differences of tlie OES parameters, 
Cfc± = (cfeo ± cj:e)/2. Columns 2-3 from the fit in) shown in 
Table II. In the last row, the proton OES has been adjusted 
by subtracting the Coulomb interaction as discussed in the 
text. Units are MeV. 



Parameter 


Cl± 


C2± 


Cfc + n 


0.88 


33. 


Ck+p 


0.93 


42. 


Ck — n 


0.06 


8. 


Cfc — p 


0.18 


17. 


^k—v 


0.10 


14. 



C. Proton fit parameters 

We now turn to the proton OES. We first note that 
the Ci+ has a slightly larger value than is found for the 
neutrons. It is tempting to interpret this as evidence 
for a stronger BCS pairing correlation energy for pro- 
tons. There is of course additional physics associated 
with the mean-field aspects of the Coulomb interaction, 
but that affects only the difference coefficient Ci-p. From 
Table III we see that the Cj_ parameters are not identical 
for neutrons and protons. The liquid-drop Coulomb en- 
ergy, proportional to Z"^ /A^^^, gives rise to contributions 

to Aop or Aep that have a very flat dependence on A. 

Parametrizing it according to model Hi) gives Coulomb 
contributions Cj+p « and ci_p « 0.08 MeV ; C2-p ~ 3. 
MeV. In the bottom row, we subtract this contribution 
from the proton parameters. One sees that c\-p is re- 



duced to a value c\_p much closer to the neutron value, 

Cl-n- 



IV. SUMMARY AND CONCLUSIONS 

The odd-even staggering of the binding energy depends 
on several effects. While a simple inverse power depen- 
dence can fit the general trends with mass number A, 
theory suggests different dependencies for the various 
terms. We have argued that three of the important con- 
tributions are single-particle energy spacing, BCS pair- 
ing, and particle-particle interactions, and that their A- 
dependence is given by A~^ for the single-particle energy 
contribution, A", i.e., constant, for BCS pairing, and A~^ 
for contributions from diagonal orbital interactions. A fit 
to the OES obtained from experimental masses indicates 
a significant improvement in the rms when both A~'^ and 
constant terms arc included. The rms values for fits with 
these two-term forms arc slightly larger than for the val- 
ues for the single inverse power form with an adjusted 
exponent. While both of these involve two free parame- 
ters the latter has no obvious theoretical basis, while the 
former does. 

The results for the constant in the fit show that an 
^-independent BCS pairing requires on the order of 
25% contribution from single-particle and time-odd in- 
teractions to the total average OES. The part left, at- 
tributable to BCS pairing, is slightly larger for protons 
than neutrons. This division implies that mean-field ef- 
fects are more important than BCS pairing in light nuclei 
{A <~ 40) and vice versa for heavy nuclei. Furthermore 
the diagonal Coulomb interaction is visible in the extract 
constants and is consistent with the liquid drop formula. 
The 1/A contributions can associated with the mecha- 
nisms of Table I, and a values for the interaction ener- 
gies can be extracted. However, this part of the analysis 
is more speculative, in view the additional assumptions 
that must be introduced. 
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